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Abstract. 

We have studied the temperature dependence of the conductivity of a silicon 
MOSFET containing sodium ions in the oxide above 20 K. We find the impurity band 
resulting from the presence of charges at the silicon-oxide interface is split into a lower 
and an upper band. We have observed activation of electrons from the upper band 
to the conduction band edge as well as from the lower to the upper band. A possible 
explanation implying the presence of Hubbard bands is given. 
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1. Introduction 

Contamination of silicon oxides by impurities has represented a long-term problem for 
Metal-Oxide-Semiconductors (MOS) production as this deeply affects the electrical 
behaviour and degrades the performance of transistors. Impurities are easily 
incorporated into silicon or the oxide throughout the fabrication of devices. The first 
source of contamination occurs during the growth of materials where the temperature 
facilitates the incorporation of impurities such as sodium, antimony, gallium and 
phosphorous into oxides and heavy metal like copper, lithium, gold into silicon. However, 
contamination also happens on processed devices that are left unprotected even in a 
cleanroom environment [1]. The silicon oxide is effectively vulnerable to various fast 
diffusers at room temperature like H + , OH~ and more specifically Na + ions. The study 
of the electronic properties of silicon transistors with oxides containing such types of 
ions is consequently of great interest from the perspective of optimising the quality of 
electronic components. Early studies started in the 70's when Fowler and Hartstein 
used a silicon MOS Field Effect Transistor (MOSFET) to probe the impurity states 
located below the conduction band in a sodium-contaminated device (2j [3]. They 
found the presence of the ions near the Si-Si02 interface at a concentration of few 
10 12 ions. cm -2 was responsible for the formation of an impurity band. At lower impurity 
concentrations, the situation is more complex as the impurity band splits into a ground 
and several excited bands [U [5] . The formation of two separate bands in sodium-doped 
MOSFETs was experimentally observed in a previous study for temperatures below 
20 K and confirmed by the analysis of the density of states at different gate voltages 
(figure 1) [6]. In the case of sodium, the ground state band is predicted to be formed by 
a single bound electron while the first excited band is formed a pair of bound electrons 
(respectively lower and upper Hubbard band). The possibility that two electrons may 
occupy the same site has been suggested by Fowler et al [7] but has remained an open 
question. This eventuality is of importance to quantum measurement based on charge 
quantum bits (qubits) (Ej. In order to investigate that possiblity, it is necessary to 
understand the electronic properties of low-doped MOS. To this end, we performed 
complemetary measurements which looked into details of the temperature dependence 
of the source-drain conductivity in the regime where electrons are activated from the 
impurity band to the conduction band edge, i.e above 20 K. We show that results are 
also consistent with an activation from the lower to the upper band as this has already 
been observed [9]. Arguments are given in favor of the presence of Hubbard bands. 

2. Experiment 

We have fabricated MOSFETs in a circular geometry (Corbino) from a (100) p-silicon 
wafer with high resistivity (10 4 n.cm). Such transistors has been widely used because 
of the ability to continuously vary the electron density and the Fermi energy by use 
of a metal gate. Its geometry eliminates leakage current paths around the contacts 
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as well as minimizes scattering with Boron acceptor impurities, especially close to 
the interface. The effective channel length and the interior diameter of the Corbino 
MOSFETs were measured to be respectively 1 and 110 /xm. A 35nm-gate oxide was 
grown at 950 °C in a dry, chlorine-free oxygen atmosphere. Contacts were realized by 
implanting phosphorous at high dose and sputtering aluminium. The contact resistivity 
was measured to be 3.5 and 2.3 fi.cm -1 respectively at nitrogen and helium temperatures 
and the sheet resistance was 6.3 and 5.9 Q.D _1 for the same temperatures. Sodium ions 
were introduced onto the oxide surface by immersing the device in a 10~ 7 N solution of 
high purity sodium chloride in deionized water. The surface of the chip was then dried 
with nitrogen gas and an aluminium gate subsequently evaporated. The application of a 
positive gate voltage (+4V at 65°C for lOmins) causes the sodium ions to drift towards 
the Si-SiC>2 interface without diffusing into silicon [TOllTT]. The application of — 4V dc 
in the same conditions removes the ions from the interface. The ions are frozen at their 
position once the device temperature becomes lower than 150 K (figurel). Standard 
low-noise lock-in techniques with an amplifier of 10 6 V/A were used to measure the 
source to drain conductivity. An ac excitation of 15 /iV and a frequency of 11 Hz were 
chosen. The dc offset of the amplifier was cut using appropriate RC filters. Finally, 
the gate voltage was controlled by a high resolution digital to analog converter and the 
temperature measured by a calibrated germanium thermometer. 

Several devices were processed identically and gave results that lead to identical 
conclusions although we noticed some variations in the relative positions and in the 
widths of impurity bands as well as in the conductivity values. We also fabricated a 
number of control devices that were not exposed to sodium contamination and were 
used for comparison. The subthreshold to saturation current ratio was about 10 3 at 
300 K and 10 6 at 0.3 K. In the case of doped devices this ratio falls down to 10 at 300 K 
and 10 3 at 0.3 K. The presence of charges close to the Si-Si02 interface is responsible 
for the finite conductivity in the accumulation region of the MOSFET. The following 
results are presented for a specific device that was chosen for its high reproducibility in 
time as well as for its high signal to noise ratio. 

3. Results and discussion 

We have measured the source-drain conductivity a for gate voltages V g between - 
2.5 and 0.4 V and for temperatures ranging from T = 0.3 to 100 K. The variation 
of conductivity with temperature was reproducible up to 110 K. We have previously 
shown that correlated hopping is present below 20 K in our device [6]. In this paper, 
we focus our analysis on the upper part of the temperature range where hopping is 
expected to disappear and be replaced by activated behaviour. However, we attempted 
to fit the curves using a general hopping model in which the exponential prefactor 
was temperature dependent and the density of state at the Fermi level was energy 
dependent. Extensive statistical analysis showed hopping does not occur above 20 K 
but that the conductivity follows an activated behaviour as described by the expression 
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Figure 1. a) Corbino MOSFET used in the experiment when the sodium ions are 
close to the Si-SiC>2 interface. The source is connected at the centre of the Corbino 
whereas the outside contact is used for the drain b) Schematic diagram of the density 
of states (DOS) for the present device, with a low energy (LHB) and a high energy 
(UHB) impurity band separated by a gap to the conduction band edge (CB). The 
dotted line shows the position of the mobility edge of the conduction band. 



a (T) = o"oT p e -€ / fcsT where <7o and e are functions of V g and p is an adjustable 
parameter. The exponential prefactor was taken to be temperature dependent in order 
to reflect the temperature dependence of the mobility of electrons [12]. Its dependence 
in temperature as well as the value for cro reflects the complex scattering mechanisms 
occuring at the Si-Si02 interface. Best fits were obtained for p — — 1 for which the 
Arrhenius plot of oT gave straight lines, indicating an activation mechanism for the 
electrons. 

Small-polarons also give rise to such a prefactor in the mobility in the adiabatic 
regime, where electrons hop without introducing further deformation [13l[Tlj. However, 
the highest temperature used in the present study is 100 K which is well below the 
energy of longitudinal optical phonons in silicon. The presence of such polarons is 
thus unlikely. However, this value for p is widely found in systems where the hopping 
transport is through a random potential or for gaussian density of states [15]. This 
is consistent with early studies on the same device proving the existence of impurities 
bands [6]. Such an activated behaviour for the conductivity has already been observed in 
MOSFETs devices [21 [16]. In our case, the situation is complex and depends on the range 
of gate voltage studied. For —0.5 V < Vg < 0.2 V, there exists two activation mechanisms, 
one at lower temperature characterized by an activation energy e\ and a second one at 
higher temperature with an energy e[. However, a single activation process with an 
energy 62 is found for — 2.5 V< Vg < —0.5 V. For comparison, two impurity bands were 
found respectively for — 2.5 V < Vg < -1.7 V (lower band) and for — 1.2 V< Vg <0.2 V 
(upper band) (Inset in figure 2). 

Our control devices containing no sodium showed no activated transport at any gate 
voltage. Below Vg = -0.3 V the conductivity of these devices remained close to zero for 
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Figure 2. Temperature dependence of the conductivity for Vg = 0.3V (A),— 0.15V 
(• ), -0.4V {<}), -0.65V (*) and -1.95V (O ) in the activated regime. The dotted 
lines defines the different regions of activation and their corresponding energies. The 
inset shows the variation of logu with V g for T=100, 20 and 0.3 K (from top to 
bottom). 
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Figure 3. a) Temperature dependence of the conductivity for an undoped device. The 
difference in threshold voltages between the undoped and the sodium-doped device is 
0.2 V. b) Variation of the conductivity of the undoped device with gate voltage at 
300 K and 0.3 K. 



the range of temperatures we investigated. Above Vg = -0.3 V and for temperatures up to 
50 or 70 K depending on the gate voltage, only correlated hopping as described by Efros 
and Shklovskii [17] was present (figure 3) and in agreement with Mason and Kravchenko 
experiments [18]. Effectively experiments were carried out in the insulating side of the 
metal-to-insulator transition where electron-electron interaction may be important in 
the absence of disorder induced by impurities at the interface. Moreover the size of the 
Coulomb gap was sufficiently wide not be screened by the temperature. This suggests 
the presence of activation energies e[ and 62 are related to the presence of impurity traps 
at the Si-Si02 interface. 



Activation mechanisms in sodium- doped Silicon MOSFETs 



6 




Figure 4. Variation of the activation energy in terms of Vg for ei (• ), e[ (y) 

and 62 (O ). The line represents the linear variation for Vg close to the threshold 
voltage. Energies are determined within an error Ae = ±0.5mcV for e < 25meV 
and Ae = ± 1 meV above. The inset shows the exponential variation of ei with gate 
voltage. 



3.1. Processes e\ and e' x 

The presence of two distinct activation energies e± and e[ suggest two different 
activation mechanisms. Unlike Fritzsche's observations [H] these mechanisms cannot 
be simultaneous because of the negative curvature of the conductivity in figure 2 at high 
temperature. Consequently, the source-drain conductivity is not written as a sum of 
parallel processes in this region. Therefore, this suggests that the mechanism responsible 
for e\ may disappear at a given temperature for the benefit of the mechanism e' v Figure 
4 shows the variation of the activation energy at different voltages. We observe that e\ 
has an exponential variation with V g . In MOSFETs, for weak accumulation, the surface 
potential energy is exponentially dependent in gate voltage. The variation of e\ (V g ) then 
simply reflects the change in the intrinsic Fermi energy and the gate capacitance when 
the gate voltage is varied. We can then assume that e± corresponds to an activation of 
electrons from the upper band to the conduction band edge [20]. The absence of the 
activation e\ for V g < —0.55 V may be interpreted as the existence of a conduction band 
threshold in the upper band separating the localized states in the band tail from a region 
of conducting states at the centre of the band. The upper band edge is therefore 36 meV 
below the conduction band edge. Close to the threshold voltage, the variation of the 
activation energy with gate voltage is linear, in the first approximation. By capacitance 
arguments, it is possible to express the value of e\ in terms of the distance d-^ + of the 
sodium ions from the interface : 

ei = *io -*^(Vg-Vt) (1) 
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Figure 5. Variation of the exponential prefactor in terms of the activation energy for 
ei (• ), e[ ( V ) and e 2 (O ). 
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Figure 6. Variation of the critical temperature T* with V g . Errors are given by 
the size of bars. The drawing shows two adjacent impurity potentials separated by 
d ~ 16 ± 1 nm with for the energy of the bound state. 



where is the threshold voltage. 

We find that the threshold voltage for conduction in the conduction band is 
~0.2V and that the shallowest localized states of the upper band reside at 
ei = 8.5meV below the conduction band edge. The thickness of our oxide being dox 
of 35 nm, the ions may then lie as close as 0.7nm from the Si-SiC>2 interface. This 
value is in good agreement with the earlier results from Di Maria [21] who obtained 
d-^ a + = 0.5 nm by measuring the oxide photocurrent at nitrogen temperature. We can 
also estimate the width of the upper impurity band tail from the energy range for which 
the activation process is present, giving a width close to 27meV. Concerning e[, the 
information obtained from the plot of the activation energy versus gate voltage is not 
sufficient to determine the nature of the corresponding mechanism and the variation of 
the pre-exponential factor with V g (figure 5) needs to be analysed. 
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The plot of c"oi( e 'i) (fig ure 5) is consistent with the Meyer- Neldel rule (MNR) 
[22] that gives an activation law for the prefactor. This behaviour is found across 
all ranges of activation energy with a slope corresponding to a typical energy Eq of 
10.5 me V whereas the relation is highly non-monotonic for a m . The observation of 
the MNR is commonly reported in single crystal, polycrystalline, amorphous, organic 
semiconductors as well as ionic crystals and glasses but more generally in inhomogeneous 
semiconductors [23] |2~H [25] . Interpretations on the origin of the MNR are various but 
the difference in the behaviour of croi( e 'i) anc ^ ^oi^i) makes our experimental results 
incompatible with a variation of the Fermi energy in temperature [26] . Multiphonon 
hopping has also been proposed [27] but this would imply the absence of the MNR 
for e[ < E which is not the case. In an ionic crystal, the MNR is explained by a 
restructuring of the lattice formed by the ions [28] and an activated conductivity due to 
an activation of the carrier density or of the mobility [29] . In our case, an analogy could 
be made except that sodium ions are unlikely to move below 100 K. We have verified 
this by performing thermal cycles from 0.3 K to 100 K that modify neither the height 
or the position of the conductivity peaks in figure 2. The mechanism responsible for e[ 
may be interpreted in the following way. Below a certain temperature T*, sodium ions 
form a disordered lattice and some electrons are localized by the impurity potentials 
near Si-SiC>2 interface at the ion site. The effect of temperature is to thermally activate 
the bound electrons for conduction with an activation energy e\ and the plot of croi(ei) 
gives information on the density of states. Above T*, the temperature is sufficiently 
high to delocalized the electrons from their sites allowing a percolation process to be 
more efficient and the scattering time to be reduced. Consequently, the conductivity is 
limited by scattering at the interface and the mobility activated with an activation e[. 
It is then possible to get e[ of the order of k^T or even smaller because the effect of high 
temperature is nomore to thermally activate carriers but to increase their mobility. This 
typically describes a process of relaxation happening at the interface like those occurring 
in melting glasses. As consequence, the slope of a' 01 (e[) may give the temperature T g 
at which the glass transition occurs [30] and the sodium ions start moving around their 
positions. We find T g ~ 122 K. The critical temperature T* could then correspond 
to the necessary energy to delocalize the electrons [31]. It is experimentally accessible 
by determining the temperature for which a change in the curvature of o~T happens in 
figure 2. Figure 6 presents the variation of T* in terms of gate voltage. As expected, the 
value of T* decreases when the electron density and the derealization are increased. 

3.2. Process ei 

The presence of e<i is explained by the finite value of the conductivity in the accumulation 
region of the MOSFET. In undoped devices, the conductivity in the same region was 
close to zero even at 100 K. Consequently, we did not find any measurable activation 
energy in the accumulation region in the reference devices. In the doped device, the 
average value for 62 is too small compared to the activation energy E\ to be related 
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to an activation to the conduction band edge and it is mostly constant over the band 
gap. It may however results from an activation to an upper state below the conduction 
band that may be inside the upper impurity band. For V g < —0.5 V, the value of 
€2 ~ 15 meV is nearly independent of the gate voltage except in the small transition 
region around -0.55 V. It is then hardly conceivable that the variation of £2 is linked 
to any bandwidth. This behaviour can possibly be explained by considering that the 
Fermi energy for V g below -0.55 V is pinned due to the difference in the number of 
states in the gap and in the band tails [32]. In that case, the Fermi energy then varies 
abruptly as V g is made more negative and gets pinned in the lower band. Also, the 
presence of a band to band activation suggests that the lower band is insulating and 
that the upper band has a conducting region. In the case of sodium ions in the silicon 
oxide, the lower band states would then be fully occupied and associated with one 
bound electron per ion whereas the upper band states would be associated with two 
electrons per ion. This situation is typically found in Hubbard bands [33J. It is then 
possible to consider that the presence of 62 is due to an electronic transition from the 
lower Hubbard band to the upper Hubbard band and that the 15 meV corresponds to 
the Hubbard gap. To strenghten this assumption, we need to compare this value to 
the theoretical estimates in the case of our device. According to the theoretical results 
of Bethe [M] on the D~ state (i.e. two electrons bound to a donor), an activation 
from the neutral state (sodium ion with a single electron) to an excited Na~ state is 
possible. Taking 11.7 for the relative permittivity of silicon as well as 0.19 m e - for the 
effective mass of electrons in silicon, we obtained €2 = 17.8 meV using Bethe's formula 
for a single D~ level. Because of the formation of a band the previous estimate for 
the activation is lowered [35] and we obtain 62 ~15meV by taking the value of the 
localization length £ ~ 23 nm measured previously on the same device [6] and a density 
of neutral donor N ~ 1.3 x 10 11 cm -2 . This value is reasonable knowing the average 
density of donor estimated from the threshold voltage shift gives Nd ~ 3.7 x 10 11 cm -2 . 
Using Nishimura's estimate, we also found the upper states in the D~ band are 3.9 meV 
below the conduction band. This value has to be compared with the 8.5 meV found for 
ei at the threshold voltage. The difference typically gives the position in energy of the 
conduction band edge relative to the bottom of the conduction band. In order to bring 
further arguments in support of the existence of a D~ band, it is necessary to give an 
estimate of the value of the on-site Coulomb energy U using the same method as Schiff 
but in two dimensions 1361 : 



e 2 

U = — (2) 

where £ is the localization length. 

The value of the localization length as well as its variation in gate voltage was 
calculated for the same device under the same experimental conditions [6]. Below 
V g = —0.5 V, £ ~ 22.5 ± 0.7 nm and the on-site repulsion energy is U ~ 5.2 meV. 
Supposing a linear relation between the gate voltage and the activation energy, it is 
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possible to convert gate voltages into energies in the gap region at a rate of 7.7meV/V 
(figure 4). The energy U then corresponds to a difference of 0.67 V in the gate voltage. 
This agrees well with the value of the existing soft gap between V g = —1.9 and 
— 1.2V in figure 2. We can also check the validity of our assumptions by calculating 
independently the value of the Hubbard gap. For sufficiently low electron density, T* 
and T g are expected to have similar values. The Hubbard gap corresponds to the energy 
necessary to put a second electron on the same site, that is E + U. We find a value 
of 15.7 meV in agreement with the value of 62- These observations are thus compatible 
with the formation of Hubbard bands and a Mott-Hubbard gap formed by Coulomb 
interactions. Finally, we would like to point out that the formation of D- states in our 
device may result from a complex interplay between inter-site Coulomb interactions, 
sodium density and disorder. The inter-site interactions are due to Coulomb interaction 
between trapped electrons in silicon. These electrons are localized in silicon at the 
potential minima created by the ions in the oxide. Disorder comes from the fluctuations 
in position of these minima in the plane but also from the height of the potential 
wells. Recent studies in GaAs/AlGaAs devices with a silicon d-doped layer showed that 
disorder could be controlled by the position of the dopant layer relatively to the 2D 
electron gas [37]. This indicates that the disorder strength in our device is at least 
partly related to the position of the ions relatively to the interface. 

4. Conclusion 

We have observed three activation mechanisms in sodium-doped silicon MOSFETs for 
20 K< T <100K. These results are consistent with the existence of two Hubbard- 
like bands, one with one electron per site and the second with two electrons per site 
|33j . The first mechanism is an activation of electrons from the upper band edge to 
the conduction band edge. The second takes over the first one when the temperature is 
sufficiently high to delocalize electrons and corresponds to the activation of the source 
drain mobility. The last mechanism has been attributed to an activation of electrons 
from a lower band to an upper band. The theoretical expectations for the position of 
the bands in energy as well as the soft gap in a case of a D~ band well agree with the 
values obtained experimentally. It is thus likely that the transport in such a localized 
system could be explained within the Hubbard model and that the observed upper and 
lower impurity bands could correspond to Hubbard bands. 
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